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Formation mechanisms of single, twin, and double hypernuclei from E~ absorption at rest 
on C are investigated with an refined microscopic transport model, that incorporates the 
recently developed Quantal Langevin treatment into Antisymmetrized Molecular Dynamics. 
The quantum fluctuations suppress the formation probability of double hyperfragments to 
around 10%, which is comparable to the experimental data, and the dynamical formation of 
twin hyperfragment can be described qualitatively. 
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§1. Introduction 



Among various nuclear fragmentation processes, the hyperfragment formation 
from the E~ absorption reaction at rest is of primary importance for understand- 
ing strangeness in nuclei. First, the S~ absorption is the most effective and the 
most direct way to produce double A nuclei. For example, all three double hyper- 
nuclear formation events in which AB^a is extracted have been discovered through 
S~ absorption on light nuclear targets. ^ " 3 ) Double hypernuclei give us valuable 
information of the low-energy YY interaction, and it might be a doorway to study 
multi-strangeness systems such as strange baryon matter, which is expected to be re- 
alized in neutron stars. 4 ) Because of this importance, further experimental searches 
for double hypernuclear formation through S~ absorption are being carried out and 
planned at BNL and KEK. 5 )" 7 ) In these experiments, the number of stopped S~ 
is expected to exceed the existing data by an order of magnitude. Therefore, the- 
oretical studies of the S~ absorption reaction and the resulting fragmentations are 
urgently needed. Second, in the KEK E176 experiment, one finds interesting frag- 
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mentation patterns, such as 12 C + £~ — > + 9 A Be, in which two single hypernuclei 
are formed. 8 **' 9 ) This kind of reaction is called twin single hypernuclear formation. 
It can be regarded as a fission of a S = —2 system or an exotic decay, which is 
unexpected in light nuclei, therefore it may suggest some special characteristics of 
strangeness in nuclei. 

In KEK E176, the estimated number of E~ absorption events at rest on light 
targets is 31.1 ± 4.8, among which one double hyperfragment formation event and 
two twin single hyperfragment formation events were observed. 10 - ) In addition, they 
found one event in which S = —2 hadronic system is formed (event type (B) in 
Table I), eight single hyperfragment formation events, and eight other events, in 
which the visible energy release is greater than 28 MeV (event type (C) in Table I). 
By using these experimental observations, the lower limits (90% confidence level) 
of S = — 2 and S < — 1 sticking probabilities were estimated as 4.8% and 47.6%, 
respectively. 10 ) In a similar way, we can estimate several lower and upper limits of 
hyperfragment formation probabilities, as shown in Table II. 

Table I. Observed stopped S~ events captured in light nucleus. In events (A), the hyperfragment 
type (double, twin, or single) is specified. In events (B), a hadronic system with S = —2 is 
formed (double or twin), and at least one hypernuclear weak decay must occur in events (C) 
(S = -1 or -2). 





Hyperfragment 


(A) 


(B) 


(C) 


S = -2 


Double 


1 


} +1 






Twin 


2 


}*■ 


S = -l 


Single 


8 







Table II. Estimated lower and upper limit (90% confidence level) by using the number of events 
shown in Table I. As for the roughly estimated values, see the text. 



Hyperfragment 




L. L. (%) U. L. (%) 


Rough 


Est. (%) 


Double 
Twin 


0.66 


} 48l0) | 47.6 10 ) 8L5 


3 - 
6- 


- 9 

- 18 


Single 


14.5 




26- 


- 73 


No Hyp. Frag. 




48.8 







These upper and lower limits are values which we must respect. However, the 
statistics is so poor that the constraints on the theory is very loose at present. 
Therefore, in this paper, we use very crude but plausible estimates of double, twin, 
and single hyperfragment formation probabilities as follows. One extreme would be 
to consider that there were no more double and twin hyperfragment formation other 
than the specified ones. Then the double, twin, and single hyperfragment formation 
probabilities can be estimated as 3% (=1/31.1), 6% (=2/31.1) and 26% (=8/31.1), 
respectively. Another extreme is to assume that all the weak decay accompanied 
by charged particles are observed, and the ratio among the double, twin and single 
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hyperfragment formation is kept. Under this assumption, the probabilities might 
be estimated as 9% (= 1/11), 18% (=2/11), and 73% (=8/11), respectively. We 
should note that the probabilities of twin and double hyper fragment formation are 
comparable. 

Various analyses of the E~ absorption reaction have been carried out. These 
studies range from statistical approaches to direct reaction theories. However, no sat- 
isfactory consistent description of double and twin hyperfragment formation is given 
yet. With statistical decay models of the double hyperon compound nucleus, n ) the 
calculated double hyperfragment formation probability is too large as compared with 
that of twin hyperfragment formation, provided that the AA interaction is attractive. 
Although the relative ratio of double- and twin-hyperfragment is improved if the AA 
interaction is assumed to be repulsive, this assumption contradicts to previously ob- 
served double hypernuclei. x > Thus the experimental data suggest the action of some 
dynamical effects that cannot be mocked up by a simple escape probability of one 
A particle after the primary elementary reaction, S~p — ► AA. On the other hand, 
Yamada and Ikeda have proposed an intuitive direct fragmentation picture which 
is based on the symmetry of the target ( 12 C) wave function: 12 )' 13 ) When proton 
s-hole and p-hole states are created in the primary elementary reaction and a highly 
excited double hypernuclear state is formed, it is fragmented mainly at the doorway 
stage to various hypernuclei by emitting some nucleons or clusters. Following this 
picture, they obtain that the excited channels of ^ 2 Be + p, ^ 2 B + n, and ^Be + d 
are produced more strongly than other channels, and the calculated S = —2 sticking 
probability is in agreement with the KEK E176 experiment. However, this model 
still underestimates the twin hypernuclei formation probability. 13 ** 

The aim of this paper is to investigate the formation mechanism of single, twin, 
and double hypernuclei from the S~ absorption at rest on 12 C by applying a re- 
cently developed microscopic transport model that augments the Antisymmetrized 
Molecular Dynamics (AMD) model 14 ) by the effect of quantum fluctuations as given 
by the Quantal Langevin model. " 

The microscopic transport approach can describe a variety of processes, from 
fast reactions, such as quasi free processes that are mainly determined by the ele- 
mentary two-body collisions, to more central reactions where bulk collective motion 
occurs and the mean-field dynamics is important. In addition, combined with sta- 
tistical decay models, the light-particle evaporation from moderately excited nuclei 
can be described. Since we want to extract the contribution of dynamical processes 
to the hyperfragment formation from the E~ absorption, the applicability to both 
dynamical and statistical processes is an important advantage of the transport ap- 
proach. 

Among various microscopic transport models, AMD and Fermionic Molecular 
Dynamics (FMD) 1S ) take account of the fermionic nature of baryons by using anti- 
symmetrized wave packets. This feature makes it possible to describe fragmentation 
processes due to cluster and nuclear shell effects, 14 )' " 20 ) which are important for 
the study of light nuclear system, including hypernuclei. 21 ) Especially two-body col- 
lision processes, which are indispensable for the description of nuclear reactions, are 
incorporated into AMD. Because of these advantages, we adopt AMD combined with 
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statistical decay model as the starting point for studying E~ absorption reactions. 

However, even if anti-symmetrized wave packets are employed, the use of product 
wave functions, whether of Gaussian form or of the Hartree-Fock type, has certain 
generic shortcomings that may be of specific importance in the present context. In 
particular, it is usually difficult to describe the dynamical formation of ground-state 
fragments in the outgoing state. This point was already noted and some remedies 
were proposed. 14 )' 23 )' 24 ) In the hyperfragment formation problem, for example, 
Nara et al. suggested the necessity of the direct productions K~ + a —> tt° + 4 H in 
the K~ absorption reaction at rest on clusterized light nuclear targets which cannot 
be described with only two-body collision terms. 24 ) 

Another possibly important problem of microscopic transport models concerns 
the inherent fluctuations associated with wave packets. When the functional space 
of single-particle wave functions is restricted to Gaussian form and the system is 
described by a single Slater determinant, the neglected degrees of freedom are ex- 
pected to cause fluctuations of the retained degrees of freedom. A treatment of this 
inherent problem has been proposed recently by Ohnishi and Randrup in the form of 
the Quantal Langevin model. 15 ) " 17 ) With the development of this model, they have 
discussed the importance of the fact that wave-packet wave functions are not eigen- 
states of the Hamiltonian operator and thus possess inherent energy fluctuations 
that affect both the statistical properties and fragmentation reactions in heavy-ion 
collisions. 

The above features may affect the fragmentation from the S~ absorption re- 
action. First, the observed twin hyper fragments are in their ground states or low 
excited states. Next, if the excitation energy is shared among all the degrees of 
freedom, each particle will have a smaller energy than the separation energy and 
then particle evaporation and fragmentation may be artificially suppressed in molec- 
ular dynamics calculations. For example, in the S~ absorption reaction, many large 
excited fragments like and 12 B are produced in the dynamical simulation of 
AMD as discussed later. This consideration imposes the importance of the quantum 
statistical features, which enhances ground-state fragments since internal degrees 
of freedom of fragments are likely to be frozen while the relative motions between 
fragments are agitated. 

Therefore, in this paper, we focus our attention on the effects of the quantum 
fluctuations. Specifically, we extend the usual AMD transport model to incorporate 
quantum energy fluctuations in the manner of the Quantal Langevin model. We call 
the resulting model AMD-QL. We will show that the application of AMD-QL to S~ 
absorption at rest on 12 C leads to total formation probabilities of various double 
hypernuclei of about 10%, when the AA interaction is attractive. 

This paper is organized as follows: We first briefly describe how we can incor- 
porate quantum fluctuations in wave packet dynamics in the manner of the Quantal 
Langevin model in Sec. 2. In Sec. 3, we incorporate this quantum fluctuation into 
AMD, and give the form and the strength of the fluctuation. Then we show our 
results of hyperfragment formation in E~ absorption at rest in Sec. 4. In Sec. 5, we 
summarize our work. 
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§2. Quantum fluctuations in wave packet dynamics 



As already mentioned in the introduction, the fluctuations are not fully incorpo- 
rated in the transport models describing the time evolution of wave packets based on 
the time-dependent variational principle (TDVP). In a statistical context, this was 
already pointed out in Refs. 15)- 17), 30). There are some claims that the statistical 
properties of the wave packet dynamics can be interpreted to be quantal rather than 
classical by re-scaling the temperature to fit the exact excitation energies and by 
analyzing the wave function itself through, for example, the single-particle energy 
spectrum 25 ^' 26 ). However, this kind of interpretation does not work in the case of 
fragment formation, since the fragments are described by the centroid parameters 
of nucleon wave packets. Thus the same authors of Refs. 25 )' 26 ) had to introduce 
additional fluctuations or correlations to describe fragmentation processes. 26 ^ " 29 ^ 

Up to now, several origins of fluctuation in wave packet dynamics have been 
proposed, such as the momentum width in a single-particle wave packet, 27 ) short- 
range correlation between nucleons, 28 ' or the lack of self-consistency between the 
wave packet and the mean field. 29 ) In this work, we do not specify the source of fluc- 
tuations, but instead require that the quantum statistical equilibrium is approached 
in the course of the time evolution following the Quantal Langevin model. 15 ** " 17 -* 
This approach has the merit that the equilibrium properties are ensured. On the 
other hand, since the system approaches the same equilibrium with any fluctuations 
within the Quantal Langevin model, we have to determine the form and strength of 
the fluctuation separately to describe dynamics. In this paper, we parametrize the 
fluctuation in a relatively simple form and seek to determine the allowed range of 
the fluctuation strength by applying the model to some phenomena with dissipation 
and/or fluctuation. 

In this section, we briefly describe the Quantal Langevin model and its appli- 
cation to dissipative nuclear collective motion, simulated within a simple soluble 
model, the Lipkin model. 22 ) 

2.1. Quantal Langevin Model 

The Quantal Langevin model is designed to ensure that the equilibrium proper- 
ties of the system are in accordance with quantum statistics. In the case of canonical 
and microcanonical ensembles, the statistical properties are governed by the parti- 
tion function and the microcanonical phase volume, respectively, 



where \ Z) represents a parametrized and normalized quantum state, and / dr \ Z)(Z\ 
resolves unity. From these expressions, we see that the probability to find a state 
\Z) is proportional to the statistical weight W(Z) = Wp{Z) or W E (Z). 




(2-1) 




(2-2) 
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To produce the desired equilibrium distribution <p(Z; t) oc W(Z) = exp(— F(Z)), 
it is possible to adopt the fluctuation-dissipation dynamics described by, for example, 
the Fokker-Planck equation, 

^),| +{ ,, W}RB .._ ? |_(,,_ ? M,^),,,,3) 

where {qi} are canonical variables satisfying dr = Ylidqi, D/Dt represents the 
time derivative along the classical path, and H denotes the classical Hamiltonian. 
The second relation is recognized as the Einstein relation that follows from the 
requirement that the equilibrium distribution be a static solution of Eq. (2-3). In 
numerical simulations, it is easier to treat the corresponding Langevin equation, 

-^T = « - {Qutyp.B. = Vi + Y^flijCj ■ (2-5) 

3 

where the fluctuation matrix g is related to the mobility tensor M through gg = M 
and ((t) is the white noise, -< (*(t) O(*0 > ~ = S(t — t'). Here we have ignored the 
diffusion-induced drift term. 31 ) The R.H.S. of Eq. (2-5), which includes both the 
average drift term and the stochastic diffusion term described with the white noise 
C, is referred to as the Quantal Langevin force and it gives the energy fluctuation 
for the system. This random force arises from the energy dispersion of each wave 
packet and thus has a purely quantal origin. Numerically Eq. (2-5) can be solved by 
the finite difference method in which white noise integrated over a short period At 
is treated as a random number following the complex normal distribution generated 
independently in each time step, multiplied by y/2At. 

In a practical application, we use the ansatz for the statistical weight W(Z) by 
applying the harmonic approximation. 16 ) For example, the statistical weight in a 
microcanonical ensemble W E (Z) is assumed to be a continuous Poisson distribution, 

W (Z) oc e -*/" {H/D)E/D - (2 . 6) 
We[Z) ^ 6 (E/D)\ ~ 6 r(E/D + l) ' 1 ' 

D = ^, al = (Z\H 2 \Z) -{Z\H\Zf , (2-7) 

where the ground-state energy is subtracted from the energy expectation value H. 
Here, D denotes a typical energy scale. Since this typical energy scale D depends 
on Z only weakly in many cases, the drift term is simplified as follows, 



where fin denotes a state dependent inverse temperature modified by quantum cor- 
rection. 16 ) With this form of the drift term, the Einstein relation becomes manifest. 
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In addition, once the energy dispersion o~% is given, it becomes feasible to solve the 
Quantal Langevin equation without further complication. 

In addition to the appearance of a random force in the equation of motion, we 
have to take the intrinsic distortion of wave packets into account before making any 
observation by using the wave packet ensemble. This point can be easily understood 
by considering the definition of statistical mean value of an observable O in canonical 
ensemble. 



< O y-. 



-I3H 



-J 



(Z\e-^/ 2 6e-^/ 2 \Z) 
{Z\e-^\Z) 



(2-9) 



Since the Boltzmann weight operator exp(— f5H) cannot be treated as a c-number, 
the distorted wave function \Z') = exp(—f3H /2)\Z) is different from the original 
evolving state, \Z). For example, the distortion operator exp(— (3H/2) emphasizes 
the energy eigencomponent with smaller energies and therefore the observed energy 
Hp = (Z'\H\Z')/(Z'\Z'} is smaller than H = (Z\H\Z)/(Z\Z}. 16 ^ 17 ^ 32 ^ 

When the total energy Eq is spec- 
ified, the wave packets are distributed 
according to the weight W Eo {Z) in the 
Quantal Langevin model, as is shown 
schematically in Fig. 1. Thus the expec- 
tation value of the Hamilton operator 
with the evolving state | Z) can fluctuate 
around the specified energy. The differ- 
ence between the evolving state \ Z) and 
the distorted state \Z') explains why 
this fluctuation of energy is allowed: 
When we make any observation, we 
have to use expectation values with the 
distorted state \Z') 




W F (Z) 



Fig. 1. Fluctuation of the energy expectation 
value. In the Quantal Langevin treatment 
with a specified energy Eq, the expectation 
value with respect to the evolving state \Z) 
fluctuates (if that wave packet has a com- 
ponent of the eigenstate of the energy Eo). 
However, any observation involves the dis- 
torted state \Z') = exp(-/3#/2)|Z). 



y/5(Eo-H)\Z), 
which is an energy eigenstate with the 
eigenvalue -Eo- Therefore, there is no 
fluctuation in the observed energy, while 



the fluctuation of the energy expecta- 
tion value with respect to the evolving state is necessary. We will discuss this point 
in the next subsection. 

2.2. Application to Lipkin Model 

In equilibrium, the Quantal Langevin model has been shown to work very well 
and give reasonable statistical properties, including the quantum statistical features 
of simple soluble systems and nuclei. 15 -* " 17 )' 30 ) In a dynamical context, the fragmen- 
tation processes in heavy-ion collisions have been studied. 17 ) However, since exact 
dynamical results are not available, the model has not been verified as a general the- 
ory to describe fluctuation-dissipation dynamics. In order to examine this aspect, 
it may be instructive to employ a schematic model, such as the Lipkin model, 22 ) 
to study the damping of collective motion, since this is one of the most familiar 
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dissipative phenomena in nuclear physics. 

The Hamiltonian operator of the Lipkin model, 



H = eK - X - V (K + K + + K_ K.) , (2-10) 



contains only the quasi-spin operators, 

^ N N 

K °=2^2 ( C +m C+m ~ cl m C-m) , K+ = c +m C-m , K- = {K+) ] . 

n=l n=l 

(2-11) 

Therefore, the collective subspace that couples to the unperturbed ground state |0) 
completely decouples from other subspaces. Within this collective subspace, the 
wave packet specified by the complex parameter z, 

\z) = exp(zK + )\0) , (2-12) 

describes a general product wave function which couples to |0). We study the damp- 
ing of collective motion through the time evolution of this wave packet by using the 
following Quantal Langevin equation for the canonical variable w = z y/N/ (1 + zz), 

2 g*v E 2 dU x dU N 

9= —' aE = ^ c to* c= oTW ( } 

When we solve this equation in z space, there naturally appears an iV dependence 
in the coefficient, g 2 z ~ g 2 a E /Nh, in the lowest order of z, which comes from the 
transformation between z and w. We have checked that the qualitative features of 
the results shown below do not depend on the detailed structure of this coefficient 
g%, such as its small-z dependence, if this N dependence is preserved. 

By using the Langevin Eq. (2-13), we have generated the ensemble of evolving 
paths. At each time, we have calculated the ensemble average of the operator K + 
with the distorted wave function, while keeping the evolving state as it is. The 
strength of the potential is chosen to keep the average potential effect constant. 
Namely, the following % parameter is chosen to be independent on N, 

x =Y-{N-l) . (2-15) 

When this parameter is less than unity, the potential energy surface has only one 
well, while it becomes a double- well shape for x > 1- As typical examples, we show 
the calculated results with \ = 1-5 an d 0.8 and iV=10, 30, and 100, and the initial 
condition z(t = 0) = 0.6. As for the fluctuation strength parameter, we show the 
results with go=0.3. However, the qualitative behavior does not change in the wide 
range of this parameter, 0.1 < go < 0.5. 

In Fig. 2, we show the time development of (ReK + ), which is related to the 
collective coordinate. For comparison we also show the exact solution and the results 
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Fig. 2. Damping of collective motion in the Lipkin model. Solid, dotted, and dashed curves show 
the exact solution and the results of TDHF and the Quantal Langevin model, respectively. The 
left (right) panels show the results with \ = V(N — l)/e = 1.5 (0.8), leading to a double (single) 
potential well. For comparison the bottom-right panel presents the case of the pure harmonic 
oscillator (one Gaussian wave packet in a harmonic oscillator potential) . 



of TDVP by using the above product type wave function. Since the above wave 
function describes a general product wave function, the results of TDVP is nothing 
but that of the time-dependent Hartree-Fock (TDHF) method. The exact solution 
shows interesting behavior. At x=1.5, the collective motion gradually damps at 
larger N, while for N = 10 the collective variable becomes larger again owing to 
the interference between some specific energy eigencomponents. In addition, we can 
see some hardening of the collective mode at smaller N. For example, one period 
becomes much smaller than the TDHF results. 

The Quantal Langevin model results exhibit a behavior similar to the exact solu- 
tion, although it is impossible to treat the interference between some specific energy 
eigencomponents which agitates the collective motion again after several periods. 
The average damping width seems to be well reproduced, and more than half of the 
hardening of this mode is explained. On the other hand, the TDHF result shows a 
pure periodic motion, since there is only one complex degree of freedom in this wave 
packet. If the potential strength parameter \ IS kept constant, there is no room for 
N to play any role in TDHF. Therefore, the results of TDHF shows the limit of 
N — > oo without any damping. 
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Within the Quantal Langevin model, the above mentioned damping and hard- 
ening can be explained as follows: Each wave packet \z) has its inherent energy 
fluctuation a E which is proportional to y/~N, and in the Quantal Langevin treat- 
ment, the energy expectation value of this evolving state can fluctuate around the 
specified energy with the order of this intrinsic energy fluctuation. Therefore, each 
wave packet will have different energy expectation value and frequencies and, as a 
result of ensemble average, the collective motion damps. The hardening of the col- 
lective mode is also explained similarly. In the case of x = l-5 there is a barrier in 
the potential energy surface at Re(z)=0 and the energy of the initial state is just 
above this barrier height. In the Quantal Langevin treatment, some of the evolving 
paths are trapped by this barrier and return to the neighborhood of the initial state 
much earlier than the classical motion. The first peak of the collective variable for 
N = 10 comes from these evolving states. These features may be related to the 
strange behavior of collective motion at around the bifurcation point. 

When N becomes larger, the importance of the above intrinsic energy fluctuation 
a E oc v^/V relative to the specified energy E oc N becomes smaller. Therefore the 
damping and the hardening of mode become smaller. 

In the semi-harmonic case x = 0-8 the situation is similar. The Quantal Langevin 
model describes the damping well, but since the potential energy surface is a single 
well and there is no trapping by the barrier, the behavior of the collective variable 
is very regular and the damping width becomes smaller. As a limit, we can consider 
pure harmonic motion. For example, we show in the right-bottom panel of Fig. 2 
the time evolution of (a') in the case of one wave packet in a harmonic oscillator 
potential. In this case, TDHF (equivalent to AMD in this case) gives the exact 
solution of the time-dependent Schrodinger equation, while the Quantal Langevin 
model still gives a small damping. However, this small damping occurs in a special 
situation and becomes visible only after many periods. Thus it is not a big problem 
since we are interested in the bulk behavior of dissipative systems such as nuclei. 

§3. Antisymmetrized Molecular Dynamics with Quantal Langevin Force 

3.1. Implementation of the Quantal Langevin Force into AMD 

In AMD, the quantum states are constructed by the following Slater determinant 
of Gaussian wave packets, 14 ) 

(3-1) 

XiU) , (3-2) 
(3-3) 

where Xi represents the spin-isospin wave function and the parameter i/j is (inversely) 
related to the variance of the Gaussian wave packet; both are assumed to remain 
constant in time. For nucleons we use = 0.16 fm~ 2 . We have included A and 
S~ particles in the framework of AMD and the width parameters for these particles 



\Z) = 



VAl detB 



det(|zi(r : 



l*(r;)> = ( ^ 



7T 



3/4 



exp 



B 
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are chosen to be vy = (my/mjv)^, where my and tun are masses of the hyperon 
and nucleon, respectively. This prescription makes it possible to factorize the center- 
of-mass wave function. The real and imaginary parts of the parameter {z{\ of the 
Gaussian wave packet (3-2), 

_ % 

Zi = y/v~idi + — ki , (3-4) 



represent the mean position d{ and the mean momentum k{, respectively. Applying 
the time-dependent variational principle to total wave function \Z), 

f d 

5 dt {Z\ih— - H\Z) = , (3-5) 
we obtain the equation of motion for the parameters {zi}, 

7£ = *-s*= - (3 ' 6) 



^ = -£C-|g, (3-7) 



J 





"J 



2 



c » = oi^- los ' MB - (3 ' 8 > 

Here, H = (Z\H\Z) is the expectation value of the total energy. 

We can construct the Quantal Langevin equation based on the AMD wave func- 
tions as before, although there are several points that require special consideration. 
The first problem is the treatment of the zero-point CM kinetic energies of frag- 
ments. 14 )' 18 )' 33 ) Since this fragment zero-point CM kinetic energy is proportional to 
the number of fragments and does not disappear even in the asymptotic region, it 
modifies the Q value of fragmentation and suppress fragmentation artificially. This 
fragment zero-point CM motion also affects the energy dispersion. For example, 
when one isolated fragment moves the energy fluctuation grows as a linear function 
of the translational kinetic energy of this fragment. In wave packet dynamics with 
product-type wave functions, where the fragment CM motion is described by a wave 
packet, these zero-point kinetic energies are indispensable. However, if we consider 
a superposition of the wave packets with respect to impact parameter and time in 
order that the incident wave function be a plane wave as in actual reactions, then 
the relative phases of this superposition are kept in time and the zero-point CM 
kinetic energies of the fragments will disappear in the asymptotic region. There- 
fore, the above effects are spurious and should be removed from the Hamiltonian as 
well as from the energy dispersion. Several prescriptions for subtracting fragment 
zero-point CM kinetic energies have already been proposed. 14 )' 18 )' 33 ) As for the cal- 
culation of the energy dispersion, we remove these effects by subtracting the local 
average velocities and forces acting on the wave packets, 

o-\ = Y^F'rCr^', , (3-9) 

ij 

F'i = Fi-^Y,fikFk , (3-10) 



Hi 



k 
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where denotes the "friendship" function which is used to subtract the fragment 
zero-point CM kinetic energies as described later in Subsec. 3.2, n« is the mass 
number of the fragment to which the i-th nucleon belongs, and here we use 7=1. 
By this subtraction, the energy fluctuations are correspondingly reduced, thereby 
ensuring that the fluctuations disappear when all the fragments have become cold. 

The second problem concerns the expression of the fluctuation matrix g. If 
we were to adopt a simple constant matrix g having only diagonal terms as in the 
Lipkin model, both of the translational motion and the intrinsic motion of isolated 
fragments close to their ground states would be affected by the fluctuation. The first 
point can be avoided by adopting a fluctuation matrix of the form, 

gij ~ Sij - 1/A F , (3-11) 

where Ap is the mass number of the fragment to which the nucleon belongs. On the 
other hand, the intrinsic (or relative) motion between nucleons are not affected when 
the same fluctuations are added. Then, we can avoid the second problem by requiring 
that all the matrix elements of g are similar within an isolated fragment near the 
ground state. In order to satisfy both requirements approximately, we employ the 
following matrix g that contains off-diagonal parts reflecting single-particle overlaps, 

(3-12) 
(3-13) 



where qi = J2k fik- I n this form, the diagonal part in Eq. (3T1) is replaced by a nar- 
row Gaussian, and the constant part for each fragment is replaced by a normalized 
friendship function having a larger width. With this form, the matrix elements of 
g become almost zero and the energy is unaffected by the fluctuation for fragments 
close to their ground states. Here we have introduced one free parameter go, which 
governs the overall strength of the fluctuation matrix. We determine this free pa- 
rameter later in Subsec. 3.4. By using this fluctuation matrix, the Quantal Langevin 
force is implemented into AMD as follows, 

Zi = l -Fi + (3nJ2 9ik9klF'i + 9ikCk ■ (3-14) 
n ki k 

In addition to introducing a stochastic force in the equation of motion, the Quan- 
tal Langevin model requires that the intrinsic distortion of wave packets be taken 
into account in order to project the state to the appropriate energy shell before 
making any observation. 16 ^' 17 )> 32 ) In the case where the energy is specified, the dis- 




tortion operator y 5{E — H) is very complicated. Therefore, we have employed the 
canonical distortion here. This can be done by using the following cooling equation 
to perform an imaginary-time evolution, 

z i = F' i , (3-15) 
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until the total energy expectation value reaches the given energy. Although this 
cooling process hardly influences the fragment configurations at the final stage of the 
reaction, it reduces the excitation energy of fragments significantly and the energy 
conservation is completely satisfied after this procedure. 

3.2. Effective Interactions 

We use the following effective interactions in this work. For the NN interaction, 

the effective interaction Volkov No.l 34 -* with Majorana parameter m = 0.575 is used. 

The zero-point kinetic energy (T CM ) is subtracted from the Hamiltonian in the same 

way as in Ref. 14). As a result, the total Hamiltonian is modified to the following 

form, „ n 

H = (Z\H\Z)-T CM , 

T CM = T A -a p T (A-N F ) , 

3 

2 N 



rru 



mrm 



E 



fij 



(3-16) 
(3-17) 

(3-18) 



fii = exp -v t 



(3-19) 



Here, To = 3h 2 v/2M is the zero-point kinetic energy of a fragment, A is the total 
mass number, and Np is the number of fragments appearing in the dynamical sim- 
ulation. The parameters involved in the T CM are determined as u t = 0.1 and a p 
= 0.772 in order to reproduce the experimental data of binding energies of normal 
nuclei as shown in Table III. 





Table III. 


Bindin 


5 enerj; 


;ies of single hyper nuclei 






Core nucleus 


3 H 


4 He 


7 Li 


«Be 


U B 




B(Core) 


cal.(MeV) 


10.63 


28.30 


39.66 


56.73 


74.51 


92.35 




exp.(MeV) 


8.48 


28.30 


39.25 


56.50 


76.21 


92.17 


Hypernucleus 


1H 


^He 


*,Li 


9 ,Be 


x iB 




B{ A A Z) 


cal.(MeV) 


12.41 


31.43 


45.32 


60.99 


87.15 






cxp.(MeV) 


10.52 


31.42 


46.05 


63.27 


87.58 




Sa 


cal.(MeV) 


1.78 


3.13 


5.66 


4.26 


12.64 






exp.(MeV) 


2.04 


3.12 


6.80 


6.77 


11.37 




Fragmentation to 






•^H+^He 


4 He+5iHe 






Q values 


cal.(MeV) 






3.26 


1.26 








cxp.(MeV) 






6.15 


3.55 







Other interactions (NA, AA) are assumed to be attractive Gaussian potentials, 

v AA = -90.12 exp(-0.935r 2 ) , (3-20) 
VAN = -43.62(0.1 - l.OPo- + 0.5P r ) exp(-0.935r 2 ) . (3-21) 

The range of the NA interaction is the same as that of two-pion exchange 21 ) and 
other parameters are chosen to fit the experimentally known hypernuclear binding 
energies of A R, ^He, A Li, 9 A Be, and 12 B within AMD wave functions as shown in 
Table III. The largest ambiguity lies in the AA interaction, since the information on 
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double hypernuclei is very scarce. We choose the parameter of the AA interaction 
to fit the AB AA of f A B(=4.9 ± 0.7 MeV) obtained in KEK E176 experiment. 3) ' 10) 
In Table III we also show the A separation energies and the Q values of some 
fragmentations that are relevant to this work. These energies are very important, 
since they largely affect the probabilities of A emission and decay to fragments in 
the dynamical evolution. The adopted effective interactions reproduce the observed 
separation energies within about 2 MeV, although the separation energies of light 
clusterized nuclei (|Li and ^Be) and their Q values of fragmentations (®Li — ► 3 H + 
^He, ^Be — > 4 He + ^He) are underestimated a little more because of the limitation 
of the AMD wave function. We will discuss the consequences of this underestimate 
later. 

3.3. Two-Body Collisions 

Two-body collisions are also included by using the physical coordinate to. 14 ) 
The following elementary collisions are included in our calculations, 

(1) N + N^N + N, (3-22) 

(2) N + A-> N + A , (3-23) 

(3) A + A-> A + A , (3-24) 

where N = n,p. The NN cross sections are taken from the experimental data. The 
cross sections involving A are calculated by using the Nijmegen model D. 35 ) The 
Nijmegen model D gives the coupling constants for YN interactions and the cross 
sections for N + A — > N + A, assuming a short-range core (r c ~ 0.5 fm) in all the 
baryon-baryon channels. The cross sections for A + A — > A + A can be obtained 
by extending the model D to SU(3) symmetry and keeping the hard-core radius as 
r c = 0.5 fm. 

3.4. Determination of the fluctuation matrix g through proton-induced reactions 

Until now, we have left one free parameter go which gives the overall strength of 
the Quantal Langevin force. Although this parameter does not affect the statistical 
properties at equilibrium, it gives the time scale of the relaxation to equilibrium. 
Thus it has some effect on the dynamics, especially in preequilibrium processes such 
as fast emission of particles or fragmentation. In this study, we determine the allowed 
range of the strength parameter go in a phenomenological fashion, by analyzing the 
fragment production cross sections of the proton-induced reaction on 12 C at E p = 
45 MeV, which has a similar excitation energy (43.5 MeV) as that of £~ absorption 
at rest on 12 C (39.5 MeV). 

The dynamical simulations are performed with AMD and AMD-QL up to the 
time 200 fm/c. The decay of the excited fragments produced during the dynamical 
stage are then treated by the multi-step binary statistical decay model denoted 
Cascade. 3? ) In this manner, we obtain the fragment mass distribution. 

First we investigate how normal AMD combined with Cascade works in describ- 
ing the proton-induced reaction on 12 C at 45 MeV. In the upper left panel of Fig 3, 
we show the results of AMD (histogram) and AMD plus Cascade (open circles) cal- 
culation for the fragment mass distribution of 12 C+p at 45 MeV. At first glance it 
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Fig. 3. Calculated fragment mass distribution in the reaction C+p at 45 MeV compared with 
the experimental data. The histograms show the calculated mass distribution at the end of the 
dynamical stage, while the lines with points show those after the subsequent statistical decay. 
Dashed and solid lines refer to AMD and AMD-QL, respectively, using go = 0.1 ~ 0.5. In the 
bottom-right panel, we compare the results of AMD and AMD-QL (go = 0.5) with the Cascade 
results starting from pure compound nucleus formation 13 N* (dotted line). The maximum 
impact parameter is taken to be 6.5 fm, which corresponds to a cross section of 1300 mb. The 
error bars indicate the experimental values. 36 - ) 



seems that the experimental data 36 ) is reasonably reproduced. However, in the pree- 
quilibrium dynamical stage, which is described by AMD, the fragmentation hardly 
occurs. One nucleon emission from 12 C, inelastic excitation of 12 C, and formation of 
the compound nucleus 13 N* are dominant, and most of the intermediate-mass frag- 
ments are formed at the statistical decay stage, as already discussed by Takemoto 
et al. 38 ) This description may not be reasonable by the following two reasons. First, 
the survival probability of the compound nucleus 13 N* at the end of dynamical sim- 
ulation seems too large. The life time of the compound nucleus 13 N* at E* = 43.5 
MeV is around r ~ 65 fm/c in Cascade, and the compound-nucleus formation cross 
section is at most 300 mb, assuming that the grazing impact parameter of around 
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3 fm. Therefore, at the time t = 200 fm/c when the dynamical calculation ends, 
around 5% of the compound nuclei 13 N* can survive, which corresponds to 15 mb by 
using the above estimate and well below the calculated value with AMD (~ 50 mb). 
The second is the lack of the preequilibrium emissions. Basically Cascade describes 
low-energy statistical decay of equilibrated excited nuclei and the angular distribu- 
tion of the ejectiles is isotropic, or at least forward-backward symmetric. However, 
at these incident energies the fragment angular distribution shows strong forward- 
backward anisotropy, as shown in the 9 Be/ 11 B + p experimental data, 39 ) caused by 
direct-like fragment production. 

Next we analyze 12 C + p at 45 MeV with AMD-QL plus Cascade and seek to 
determine the allowed range of the parameter go by comparing the calculated frag- 
ment mass distribution with data. In Fig. 3, we show the results of the AMD-QL 
(plus Cascade) calculation for the fragment mass distribution of 12 C + p at 45 MeV 
with 5o=0.1 ~ 0.5. After the statistical decays we cannot see significant differences 
among the results of the AMD and AMD-QL calculations, and all of them well repro- 
duce the data. However, the mass distribution at the dynamical stage (histograms 
in Fig. 3) changes drastically when the value of go is changed. As the fluctuation 
strength is increased from go=0 (AMD) to go = 0.5, the compound nuclear survival 
probability decreases and dynamical fragmentation (mainly to a+ 8 Be and 3a) oc- 
curs more frequently. We should note that AMD-QL can describe fragmentation 
to intermediate-mass fragments at low excitation before statistical decay by the ef- 
fects of fluctuations around classical paths given by Quantal Langevin force. So, the 
mass distributions before and after Cascade become closer, and direct production 
of ground-state fragments is partly described with AMD-QL. On the basis of this 
analysis of proton-induced reactions, we find the range go < 0.5 for the fluctuation 
strength parameter is suitable for reproducing the fragment mass distribution after 
Cascade. In addition, from the consideration of the life-time of compound nuclei and 
the necessity of preequilibrium fragment emission, the values around go ~ 0.5 seems 
preferable. 

It is also possible to describe the reaction 12 C + p at 55 MeV with the same 
value of the strength parameter go- Detailed analysis of proton-induced and heavy- 
ion reactions with AMD-QL will be reported elsewhere. 
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§4. Results for the £~ absorption reaction 

We now apply the AMD-QL simulation, and the subsequent Cascade process, 
to the absorption of S~ at rest on the nucleus 12 C. 

4.1. Initial Condition 

The initial wave function of the S~ hyperon is calculated by assuming the in- 
teraction between E~ and 12 C to be the following Woods-Saxon potential and the 
Coulomb potential, 

1 + U C oui(r) , (4-1) 

where R = r^A 1 ^ denotes the nuclear size, and Vq is the potential depth. Here, we 
ignore the effects of the imaginary part of Vq on the wave function. The size and the 
diffuseness parameters are taken as ro = 1.14 fm and a = 0.65 fm, respectively. 

Once the S~ wave function is known, the absorption point of S~ is calculated 
by the density overlap between E~ and the protons in 12 C, 

^oc Pp (r) te -(r) . (4-2) 

It is noteworthy that this density overlap does not depend on the radial quantum 
number n for higher I if it is normalized, since it is sensitive only to the behavior 
at the nuclear surface. We assume that S~ always reacts with a proton in 12 C and 
two A hyperons are produced at the points where the S~ and proton existed. We 
neglect the elementary reactions involving the U hyperon, since those channles are 
closed. This elementary reaction is described by ordinary collision term based on 
the canonical variable w. 14 ) The spins of the produced A hyperons are determined 
statistically. The energy of the S~ absorption on 12 C in the experiment in which 
twin hypernuclei was discovered 8 ^' 9 -* was well reproduced by a shallow potential 
(Vq = — 16 MeV) and two events of twin hypernuclear formation, 12 C + S~ — > 9 Be 
+ 4 H, are interpreted as the reactions in which S~ is absorbed from a 2p state. 40 - ) 
In this study, we assume that 5~ is absorbed from a (1 or 2)p state and Vq = — 16 
MeV. 

4.2. Hyperfragment Mass Distribution 

We have made simulation calculations by using AMD and AMD-QL (go=OA, 
0.3, 0.4, and 0.5) up to t = 200 fm/c. A few thousand events were generated for 
each value of the strength parameter g . In the case of AMD-QL, we have taken the 
effects of intrinsic distortion of wave packets into account at the end of the reaction 
by the imaginary-time evolution given in Eq. (3-15). In this procedure, we have 
ignored those events in which the total energy expectation value has not reached the 
specified energy Eq. In the statistical decay stage, the fragment mass and isotope 
distributions are determined essentially by the Q values of the decay, the Coulomb 
barrier, and the spin-degeneracy factor. Therefore, we have used the experimentally 
known hypernuclear binding energies, if available. If not, the separation energies 



U(r) = Vq 



1 + exp( 



r- R s 
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calculated by using the AMD wave functions are adopted. On the other hand, the 
fragment excitation energies used as the inputs for the statistical decay are calculated 
by subtracting the calculated ground-state energies in order to avoid the spurious 
effects arising from the difference between the calculated binding energies in the 
model space and the experimental values. 

In Fig. 4, we show the calculated single and double hyperfragment mass distri- 
bution in E~ absorption reaction at rest on 12 C with AMD and AMD-QL combined 
with Cascade. 




2 4 6 8 10 12 2 4 6 8 10 12 
Fragment Mass Fragment Mass 



Fig. 4. Calculated single and double hypernuclear mass distributions in E~ absorption reactions 
at rest on 12 C. Histograms in the upper, middle, and lower panels show calculated results with 
AMD, AMD-QL (g = 0.3), and AMD-QL (g = 0.5) after the dynamical stage, respectively. 
Open circles (upper panel), squares (middle), and solid circles (lower) show the results with 
AMD, AMD-QL (g = 0.3), and AMD-QL (g = 0.5) followed by Cascade, respectively. For 
comparison, the upper panel also includes the results for AMD-QL (g a — 0.5) followed by 
Cascade. 
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For the single-hyperfragment mass distribution, the situation is similar to the 
case of the 12 C+p reaction: With AMD, we only obtain formation of the 12 B* 
compound nucleus by one A evaporation in the dynamical stage. Almost all the 
single-hyperfragments are formed through the statistical decay of 12 B and *^B . 
On the other hand, with AMD-QL, fragmentations occur more frequently, mainly to 
^B, ^Be, ®Li, and ^He, which have relatively large binding energies and have stable 
fragmentation partners after one A emission. As the fluctuation strength go is in- 
creased, dynamical fragmentation enhances. In addition, we can see some signals of 
the dynamical effect in the single-hyperfragment mass distribution in AMD-QL fol- 
lowed by Cascade. For instance, A=10 ( 1 °Be and 1 °B) fragments are created mainly 
through one nucleon (or A) evaporation from dynamically produced ^4=11 single (or 
double) hyperfragments in AMD-QL followed by Cascade, but these processes are 
not possible in AMD plus Cascade. 

For the double hyperfragment mass distribution, the effects of quantal fluctua- 
tions of AMD-QL are more evident. For instance, the production rates of ^4=11 and 
12 double hyperfragments calculated with AMD-QL followed by Cascade are much 
smaller than those of AMD plus Cascade. This difference comes from the compound 
nucleus ^B* formation probability at the end of the dynamical stage. (Although the 
above difference looks small on a logarithmic scale, it is significant on a linear scale, 
see Fig. 6.) This formation probability is about 30% in AMD, while it decreases to 
around 10% for go= 0.5 in AMD-QL. In other words, the A emission probability is 
enhanced in AMD-QL. 

The fact that the formation of long-lived double-hyperon compound nuclei occurs 
so frequently can be understood as follows. In the initial stage of the S~ absorption 
reaction, a proton in 12 C is converted to a A and the residual nucleus becomes n B*. 
The separation energy of this proton is around 16 MeV, and 11 B* is formed with 
about 6 MeV excitation in AMD. As a result, a large part of the released energy in the 
elementary reaction of S~p — > AA (28.3 MeV) is expended and the energy available 
to the 2A system is very small (~ 3 MeV/particle). In AMD, the quantal fluctuations 
associated with wave packets are neglected and the motion of A is limited by the 
given small energy. Therefore, the formation of long-lived double-hyperon compound 
nuclei occurs more frequently in AMD calculation. 

In order to show the mechanism of A evaporation in the dynamical stage of 
AMD-QL, we analyze the dynamical evolution of the single-particle energy of each 
of the two produced A hyperons. Figure 5 displays a typical time development of 
two A hyperons in a projection showing the single-particle energy and the distance 
from the center of the nucleus. The left (right) panel is for AMD-QL (AMD) and the 
initial projections are shown by solid dots. Since in AMD-QL the Quantal Langevin 
force continually kicks the particles, one of the hyperons is occasionally emitted 
towards the continuum region, leading to its evaporation from the highly excited 
^B* compound nucleus. In AMD-QL, this type of A evaporation process dominates 
and the double hyperfragment formation is suppressed. By contrast, with the AMD 
simulation the hyperons are steadily losing their energy, leading to their eventual 
absorption into the compound nucleus. In this case two hyperons lose most of their 
single-particle energy as a result of collisions with the nucleons. In AMD this type 
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of energy loss occurs frequently and double-hyperon compound nucleus formation 
therefore becomes more frequent. 




012345 012345 
Distance From the Center [fm] Distance From the Center [fm] 

Fig. 5. Typical time development of the single-particle energies of the two A hyperons, as obtained 
with AMD-QL (left) and AMD (right). The abscissa is the distance from the center-of-mass of 
the system and the ordinate is the single-particle energy. The solid circles indicate the initial 
position of the hyperons at the time of the S~p conversion to AA. The single-particle potential 
energy of A with respect to ^ 2 B is shown by the dotted curve. In the AMD-QL simulation, 
the hyperons are kicked by the random force and one is emitted (dashed trajectory in the left 
panel). On the other hand, in AMD simulation, even the hyperon that initially had a positive 
single-particle energy is absorbed into the hyperon compound nucleus due to the energy lost by 
the collision with a resident nucleon (solid trajectory in the right panel). 



It is interesting to compare the 12 C + S~ reaction, in which two A hyperons 
are produced, with the 12 C+p reaction in which all particles are nucleons. In the 
latter case we cannot tell from the mass distribution alone whether the particle 
was emitted in the dynamical stage or in the statistical stage because the incident 
(leading) particle and the residents of the target nucleus are identical. Therefore, 
the mass distribution does not depend on whether the ejectile is the leading particle 
or arises from the target. 

On the other hand, in the S~ absorption reaction, the leading two hyperons 
are different from nucleons and we can tell whether the ejectile is one of those or 
not. During the dynamical (or preequilibrium) stage, the probability for leading 
particle A to escape is larger than that for the resident nucleons. However, there 
is no difference once equilibrium has been achieved. Therefore, we may be able to 
clarify the contribution of the dynamical fragmentation and statistical one from the 
A emission probability of the S~ absorption reaction. 

4.3. Event Type Analysis 

Table IV shows the dependence on the fluctuation strength parameter go of the 
relative abundance of different event types; double-, single-, and twin-hyperfragment 
formation, as well as events with no hyperfragment formation. The AMD results 
correspond to go=0. As described in the previous subsection, quantum fluctuations 
enhance preequilibrium emission of one A particle, while suppressing formation of 
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the double- hyperon compound nucleus (second column in Table IV). This fea- 
ture also applies to emission of other particles. In the dynamical stage, the double 
hyperfragment formation probability monotonically decreases and two-yl emission 
events (no hyperfragment formation) grows, as go increases. For single hyperfrag- 
ment formation, it first grows following enhanced one A emission, but it saturates 
and decreases again when two A emission becomes visible. 

Table IV. Dependence on the fluctuation strength go of the event type abundance (%) for double, 
single, twin, and no hyperfragment formation. Pure statistical decay results with Cascade and 
the statistical multifragmentation model are also shown for comparison. In the results of the 
statistical multifragmentation model, the sum of single and no hyperfragment formation rates 
is shown with superscript a. 





Dynamical(%) 








+Cascade(%) 






go 


13 rj* 

aa£> 


Double 


Single 


Twin 


No Hyp. 


Double 


Single 


Twin 


No Hyp. 


0.0 (AMD) 


30.3 


30.3 


66.3 


0.0 


3.4 


18.1 


39.7 


0.9 


41.3 


0.1 


29.6 


30.1 


66.1 


0.0 


3.8 


18.1 


38.4 


0.8 


42.7 


0.3 


21.0 


24.1 


71.3 


0.1 


4.5 


15.1 


38.0 


0.7 


46.1 


0.4 


16.4 


20.4 


72.0 


0.3 


7.3 


12.9 


42.2 


0.9 


44.0 


0.5 


9.2 


16.3 


70.9 


0.9 


11.9 


11.4 


44.0 


1.2 


43.4 


Cascade 


100.0 


100.0 








58.8 


24.5 


3.1 


13.6 


Compound 11 ' 


100.0 


100.0 








66 


20 a 


14 


a 



As the fluctuation strength go is increased, the probabilities after the statistical 
decay stage generally become closer to those after the dynamical stage. This feature 
is understood as arising from the quantum statistical nature: In AMD-QL with a 
sufficient fluctuation strength, a large part of dynamically produced light fragments 
have excitation energies that are sufficiently small to ensure survival during the 
statistical decay stage. This is because the lower energy components in the wave 
packet are emphasized by quantum statistics, as realized in the calculation of the 
intrinsic distortion of wave packets. Thus the difference between total production 
rates of hypernuclei at the end of the AMD-QL simulation and after the subsequent 
statistical decay becomes smaller, while the Cascade after-burner is still important 
in describing the decay of heavier fragments, such as ^B, \^B, ^Be and 12 A B. 

The combined results of AMD and AMD-QL with Cascade show that all of 
these probabilities are in the allowed range between the lower and upper limits 
experimentally estimated at 90% confidence level, shown in Table II, and the single 
hyperfragment formation probability is in the roughly estimated range (26-73%). 
These things apply at any go value. In addition, with a larger fluctuation strength, 
the double (twin) hyperfragment formation probability decreases (increases) and 
becomes closer to the rough estimate 3-9% (6-18%). 

In Fig. 6, we show the calculated double- and twin-hypernuclear formation prob- 
abilities as a function of the fluctuation strength go- At small values of go, double 
hypernuclei are calculated to be formed with a large probability of around 20%. 
When we incorporate quantum fluctuations, this probability decreases to around 
10% at 5o=0.5 due to the enhanced emission of A in the dynamical stage. On the 
other hand, the twin formation probability first becomes smaller at small go values 
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Fig. 6. Double- and twin-hypernuclear formation probabilities as functions of the fluctuation 
strength parameter g . Dashed and solid lines connect the results before and after statisti- 
cal decays, respectively. The dotted curve in the left panel shows the probability of forming 
the double- hyperon compound nucleus (^B ), and the dotted curve in the right panel denotes 
the probability of twin hyperfragment formation purely from the statistical decay of the excited 
double hypernucleus. The experimental estimates are shown by thin solid lines. 



and then grows again for g$ > 0.3. Within the model calculation made here, this 
behavior is governed by the competition between two different formation mechanism 
of twin hypernuclei; the statistical formation from the double-hyperon compound nu- 
cleus and dynamical fragmentation. At small go values, the former dominates twin 
events and no dynamical formation of twin hypernuclei can be seen. The probability 
of the latter grows as a linear or quadratic function of go, and becomes dominant at 
large quantum fluctuations, although the twin hyperfragment formation probability 
is still significantly smaller than the rough estimate. 

In Fig. 7 we show the density evolu- 
tion of a twin hyperfragment (^Li+^He) 
formation event obtained with the 
AMD-QL dynamical simulation. In 
the initial state, a highly excited ^B* 
compound-like state is formed. 

Around 80 fm/c the compound sys- 
tem starts to fission, leading to the for- 
mation of the hypernuclei ^Li and ^He 
at very small excitation. This twin 
fragmentation occurs over a relatively 
long time scale between direct and sta- 
tistical decay processes. In addition 
to twin hyper fragments, frequent light- 
fragment emission (such as 4 He, 3 H, and 
7 Li) appear in the dynamical stage of 
AMD-QL. With quantum fluctuations 
included we can obtain this variety of 




Fig. 7. Time evolution of the matter density in 
a twin hyperfragment formation event, 12 C 
+ S~ — > ® Li + ^He, calculated with AMD- 
QL. Crosses and triangles indicate the A 
and proton positions, respectively. Neu- 
trons are omitted for simplicity. 
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fragmentations during the dynamical stage, before the statistical decay. 
4.4. Exclusive Channels 

Table V shows the calculated rates for various channels with AMD and AMD-QL 
(do = 0.5), both followed by Cascade. We also show the calculated production rates 
obtained by applying direct reaction theory to the absorption from the 2p atomic 
state of S~ with Vq = —16 MeV. 13 ) Although they include the production rates to 
the particle unbound excited states, it is still valuable to compare the trend. 

As can be seen in Fig. 4, double hypernuclei are mainly produced at large masses 
(A=W, 11, and 12), because of their larger 2A binding energies. This feature partly 
explains why it is difficult to detect double hypernuclei, since the non-mesic decay 
rate is expected to be dominant in these heavier double hypernuclei. One exception 
is ^He (lampha), whose production rate reaches around 1.5%. While we can see 
some visible differences between the results of AMD plus Cascade and AMD-QL plus 
Cascade, for the reason described in the previous subsection, it is natural that the 
channels with large Q values are populated abundantly, since a large part of these 
double hypernuclei are produced through the statistical decay of excited double 
hypernuclei, ^B*, ^B*, f A Be*. 

As for the single hypernuclei, those with mass number A=ll dominate. Once 
a A is emitted, the residual nucleus 12 A B, which has the largest probability in the 
dynamical stage, will be excited less than 25 MeV, depending on the energy of 
the emitted A. These excitation energies are comparable to the nucleon and A 
separation energies and therefore a single emission is enough to bring the residue 
into the particle stable regime. The abundant single hypernuclei which are produced 
from 1 ^B, 1 ^B, 1 ^Be, ^Be, ®Li, and ^He, have their stable decay partners with small 
separation energies from ^B, S n = 12.6 MeV, S p = 14.1 MeV (AMD calc), 5 t =15.8 
MeV, S a = 13.2 MeV, and S( A He)=16.9 MeV, respectively. In AMD-QL, ^Be has 
also a relatively large probability of 1.52%, which is produced through one-nucleon 
emission from dynamically produced U A B and 1 ^Be, and deuteron emission from 12 A B. 
The dominance of heavier single hypernuclei, which mainly decay in non-mesic ways, 
explains again why we frequently see those events with large energy release without 
pion (8 events out of an average of 31.1). 

In twin hypernuclear formation channels, the difference between AMD and 
AMD-QL appears most clearly. While all the twin hypernuclear formation is de- 
scribed by the statistical decay of A 3 A B in AMD, it is mainly described in the dy- 
namical stage in AMD-QL with a sufficient fluctuation strength. Especially, in the 
channels ^Li+^He, 2^He+t, and ^He+^H+a, most of the dynamically produced 
fragments are bound and therefore survive the statistical decay stage. Another visi- 
ble channel in the dynamical stage is ^Be+^,H, which is the experimentally observed 
one. However, in the ^Be+^H channel, dynamically produced ^Be hyperfragments 
are usually in their particle-unstable states which decay into ^He+a. One reason for 
this instability is the underestimate of the separation energy to the ^He+a channel 
within the present effective interaction and AMD wave function. Another reason 
might be in the extended density profile of ^Be. The fluctuation matrix with the 
form of Eq. (3-12) is selected so as not to affect the intrinsic motion of well-isolated 
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Table V. Calculated production rates with AMD and AMD-QL (go = 0.5), followed by Cascade, 
for double, single, and twin hypernuclear formation. The results of direct reaction theory 13 - ) are 
also shown for comparison. For the value of + 2a in direct reaction theory we take the rate 
of the AA R+ 8 Be channel. 



Channel AMD " AMD-QL(p = 0.5) Direct 13) 

Dynamical(%) +Cascade(%) Dynamical(%) +Cascade(%) 



Double Total 


30.3 


18.07 


16.3 


11.39 




13 T3 


30.3 




9.2 







^B+n 





2.00 


2.4 


0.61 


3.96 


WBe+p 


0.03 


0.96 


2.1 


0.34 


2.43 


"^B+2n 





1.57 


0.4 


1.00 




" A Be+d 





0.26 


0.3 


0.09 


0.59 


"ziBe+pn 





5.27 


1.0 


2.88 




\°a Be+t 





0.12 


0.4 


0.03 


0.14 


^Be+(dn,pnn) 





2.91 


0.2 


3.17 




A ®Be+tn 








0.1 


0.10 




AA 9 Li+a 





0.16 


0.2 


0.05 


0.05 


AA Li+cm 





0.51 





0.19 




yi^He+ap 





0.36 





0.24 




AA He +( 6 Li,ad,apn) 
^He+ 7 Li 





0.79 





0.49 







0.23 





0.13 


0.05 


^He+ 6 Li+n 





0.60 





0.37 




yi^He+a+(t,dn,pnn) 





1.85 


0.1 


1.39 




AA B+2a 





0.29 





0.09 


0.11" 


Other 





0.19 




0.23 




Single Total 


66.3 


39.68 


70.9 


43.96 






66.3 


1.42 


53.0 


3.24 


2.07 


u A B+A+n 


_ 


17.79 


7.8 


20.90 




^Be+A+p 


— 


5.39 


3.7 


5.31 




10 A B+A+2n 




0.01 




0.01 




1 °Be+/l+(d,pn) 




0.79 


0.5 


1.52 




A Be+A+t 




4.20 


0.7 


2.73 




A Be+A+dn 




0.03 




0.02 




*Li+yl+a 




3.67 


2.9 


4.96 




7 A Li+A+cm 






0.1 


0.10 




A tte+A+ 7 U 




6.26 


1.1 


3.99 




5 A Re+A+at 




0.03 


0.7 


0.72 




4 A U+A+(2a, 8 Be) 




0.09 


0.4 


0.48 




Twin Total 


0.0 


0.92 


0.9 


1.22 




^Be+ 3 H 




0.08 




0.01 




'Be+lH 




0.16 


0.3 


0.04 


0.22 


^Be+ 3 H+n 




0.01 




0.01 




»Li+ 5 A He 




0.06 


0.3 


0.33 


0.21 


7 A U+ A Ue 




0.05 




0.01 




^Li+^He+n 




0.08 




0.07 




2^He+t 




0.36 


0.2 


0.34 




iHe+lH+a 




0.12 


0.1 


0.41 
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compact fragments. However, ^Be has an extended aaA structure, and so the in- 
trinsic motion will be affected especially when fragments are separating. 

The total probability of twin hypernuclear formation shown here is 0.9 - 1.2%, 
depending on the fluctuation strength. When we limit the initial system to be 
the double-hyperon compound nucleus , the Cascade probability resulting in 
twin hyperfragments is only around Pr w i n (^B*)= 3%. On the other hand, this 
probability amounts to PrwinQ^B ) =16%, within the hyperon compound nucleus 
picture, n ) which exceeds the Cascade result by a factor of five. We have found 
that it is suppressed to be around frwinQ^B ) = 10% and 7% within a statistical 
multifragmentation model similar to that used in Ref. 11), when we take account 
of excited levels of daughter fragments with mass number A > 5 up to either E* = 
E'thr + 2(^4 — 4) MeV or infinity, respectively. Therefore, a factor of two may be 
explained by a different treatment of the excited fragment states at the statistical 
decay stage. The remaining factor 2-3 may be due to differences in decay scheme. 

§5. Summary and Discussion 

We have investigated single, twin, and double hypernuclear formation from S~ 
absorption at rest on 12 C by using Antisymmetrized Molecular Dynamics 14 ) aug- 
mented by the Quantal Langevin force 15 ) " 1? ) (AMD-QL) and followed by the multi- 
step binary statistical decay model (Cascade). 37 ) 

The Quantal Langevin treatment ensures that the dynamical evolution of the 
system leads towards quantum statistical equilibrium. Such dynamics can be de- 
scribed by the time evolution of a distribution function that satisfies a Fokker-Planck 
equation, and the equivalent time evolution of the complex wave-packet parameter 
{zi} is governed by a Langevin-type equation of motion. This model has been shown 
to work well in a statistical context; it describes the statistical equilibrium properties 
of simple soluble models and finite nuclei well, 15 ) " 17 ) and it has been used to study 
equilibrium properties of argon fluids. 32 - ) Furthermore, it has met with some success 
in describing the fragment yields in Au+Au collisions at higher energies (100 - 400 
MeV/A) where the number of particles are large and the statistical equilibrium may 
be achieved. 17 ) Therefore, the model may be expected to work well also in dynamical 
contexts where pre- or non-equilibrium aspects are important. 

In order to show the validity of the Quantal Langevin model in a dynamical 
context, we first applied it to the damping of collective motion within a simple 
soluble Lipkin model. It was shown that the damping process is well described 
within the Quantal Langevin model, including its dependence on the quasi-spin 
degeneracy N. Although it is difficult to treat the interference between specific 
energy eigencomponents that agitate the collective motion again after several periods, 
the early evolution coincides well with that of the exact solution. 

On the basis of the above success in the Lipkin model, we then incorporated 
the Quantal Langevin force into the Antisymmetrized Molecular Dynamics (AMD). 
In the resulting model, referred to as AMD-QL, the usual AMD equation of motion 
for {z^ is augmented by the inclusion of the quantum fluctuations, leading to a 
Langevin-type equation. In order to achieve reliability in the practical applications, 
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we have employed a relatively simple parametrized matrix g and determined the 
range of the overall fluctuation strength parameter, go, by analyzing the fragment 
mass distribution of 12 C+p reaction at 45 MeV, which has an excitation similar to 
that of the reaction under consideration, E~ absorption at rest on 12 C. This analysis 
shows that although the fragment mass distribution of AMD-QL resembles that of 
AMD after statistical decays, it is significantly different in the dynamical stage. In 
AMD, inelastic excitation of 12 C, compound nuclear formation, and one-nucleon 
emission exhaust almost all the events and no fragmentation can be seen. On the 
other hand, in AMD-QL with a sufficiently large fluctuation strength go, various 
fragmentations occur dynamically and the mass spectrum in the dynamical stage is 
closer to the final post-Cascade results as go is increased. Considering the lifetime of 
the compound nucleus 13 N* and the predominance of preequilibrium fragmentation 
(which is evident from the forward-peaked angular distribution of the fragments, as 
is apparent in the data for the 9 Be/ 11 B + p reactions at 45 MeV, 39 ) for example), it 
appears that this new framework is superior to AMD, since the incorporation of the 
quantum fluctuation of the energy leads to a significantly improved description of 
particle evaporation and fragmentation in dynamical process. AMD-QL can describe 
the fragmentation of the moderately excited nuclei dynamically. This feature is 
expected to be especially important for the E~ absorption reaction, because the total 
excitation energy in this reaction is very small and a large part of stochastic two-body 
collisions are blocked. The quantum fluctuations should then be the primary source 
in the system for generating the fluctuations leading towards the various channels in 
the final state. 

We have employed AMD-QL for studying the dynamical mechanisms for A emis- 
sion and twin hyperfragment formation in the S~ absorption reaction. While AMD 
(plus Cascade) gives large production probabilities for double hypernuclei (18%) and 
cannot describe twin hyperfragment formation in the dynamical stage, AMD-QL 
(plus Cascade) suppresses the formation of double hypernuclei (~10%) and achieves 
a qualitative description of twin hyperfragment formation in the dynamical stage 
(~0.9% in dynamical stage, and around 1.2% after Cascade). The calculated proba- 
bilities of double, single, and twin hyperfragment formation are in the region allowed 
by the experimental estimate at 90% confidence level shown in Table II. Moreover, 
the calculated single hyperfragment formation probability is within rough estimate, 
26-73%. In addition, the double hyperfragment formation probability becomes very 
close to the rough estimate, 3-9%, when the fluctuation strength is sufficiently large. 
This is mainly due to the enhanced A emissions. In AMD, even when a A has suffi- 
cient energy to escape the compound nucleus, it easily loses its energy by two-body 
collisions before getting out. Therefore, in AMD a A can only escape from the nu- 
cleus when it is created on the nuclear surface and its initial momentum is directed 
outwards. By contrast, in AMD-QL the quantum fluctuation continually kicks all 
the particles and a A can then be emitted even if it is initially bound. Through 
these analyses, we conclude that enhanced A emission in the preequilibrium stage is 
essential in describing the bulk dynamics of S~ absorption at rest and the quantum 
fluctuation given by the form of Eq. (3-12) describes this enhancement well with the 
overall quantum fluctuation strength parameter go being around 0.5. 
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Fig. 8. Various model analyses of relative abundance of the various event types: double, twin, 
S — — 2, single, S < —1, and no hyperfragment formation. Solid triangles and solid squares 
show the results calculated with AMD and AMD-QL (go = 0.5), and the error bars in AMD-QL 
arise from the fluctuation strength parameter, from go =0 (AMD) to go=0.5. The calculated 
results based on the double-hyperon compound nucleus picture of (H.C.) with attractive 

A- A interaction 11 - 1 and the direct reaction picture for S~ absorption from 2p atomic orbit 13 ^ 
are shown by open squares and open circles, respectively. The arrows from H.C. and direct 
reaction results show the model uncertainty according to the A emission in the preequilibrium 
stage and multistep processes other than direct break-up, respectively. Experimentally allowed 
ranges are closed by solid lines (90% confidence level) and dashed lines (rough estimate). 



Although the calculated twin hyperfragment formation probabilities are larger 
than the experimental lower limit (0.66%), it is still much below the rough estimate, 
6-18%. Furthermore, both of the two twin hypernuclear formation events are in the 
channel ^Be+^H and it is known that these twin hypernucleus are formed after S~ is 
absorbed in 12 C. 8 )' 9 ) Therefore the probability for finding the twin hyperfragments 
9 A Be+\R would reach around 13% (=2 events/ (31.1 events (C,N,0) x 0.48 (C))) 
when the S~ is absorbed in 12 C. The calculated probabilities for this channel are 
only 0.16% and 0.04%, with AMD and AMD-QL (g =0.5), followed by Cascade, 
respectively. Even if we sum up the probabilities in the channels ^He+^H+a (0.12% 
and 0.41%) and 4 A K+A+(2a, 8 Be) (0.09% and 0.48%) to compensate the underesti- 
mate of the separation energy in A Be, the total probabilities, 0.37% and 0.93%, are 
still smaller than the above rough estimate by an order of magnitude. 

This underestimate of twin hyperfragment formation is not a specific problem in 
transport models, and it was already pointed out in previous works. As pointed out 
by Yamamoto, Sano, and Wakai, n ) the statistical multifragmentation decay model 
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of the double-hyperon compound nucleus also fails to reproduce the experimental 
fact that the twin hyperfragment formation is more frequent than that of double, 
assuming that the A- A interaction is attractive. Yamada and Ikeda 13 ^ also showed 
that the channels of twin hypernuclei gave very small spectroscopic factors, and the 
calculated formation probability of twin hypernuclei was very small. 

Figure 8, summarizes the results of various model analyses. Since the emission 
probability of one A at the preequilibrium stage is not taken into account in Ref. 11), 
we assume that the calculated values give the upper limit of the model. On the other 
hand, the probabilities shown in Ref. 13) for double, twin, and S = — 2 in the direct 
reaction picture are considered to be lower limits, since the dissipation of the doorway 
state into compound states and formation from these compound states are neglected. 
We can see immediately that all the analyses underestimate the twin hyperfragment 
formation probability relative to that for double hyperfragment formation. 

Therefore, if the twin hypernuclear formation probability, especially in the chan- 
nel ^Be+^H, is really high as suggested, there must be some production mechanism 
which is missed in all the theoretical works, including this work. Another possibility 
is that the twin hypernuclear formation probability is very close to or below the lower 
limit at 90% confidence level estimated by using the number of observed events. In 
order to distinguish these two possibilities, higher statistics experiments are required. 
In addition to the higher statistics, one of the key challenges is to narrow down the 
upper limit in each type of events. This requires resolving the unspecified events, 
such as those events of type (C) in Table I. 

We have noted the shortcomings of the standard microscopic AMD model with 
which it is hard to describe production of fragments near their ground states in low- 
energy preequilibrium dynamical processes, such as twin hypernuclear formation. 
As the inclusion of quantum fluctuations leads to a significant improvement, the 
extended model may yield a better description of other fission processes as well. 

It is obviously important to broaden the confrontation of the theory with ex- 
periment through application to other reactions for which suitable data is available. 
Moreover, it should be kept in mind that there are still ambiguities in AMD-QL. 
The key quantity concerning fluctuations, the matrix g, is determined phenomeno- 
logically. Therefore, in the future, it would be of interest to determine g from more 
fundamental considerations. In this context, we draw attention to a recently devel- 
oped extended AMD (called AMD-V) that incorporates a diffusion of single-particle 
wave functions. 29 ) Since this description also introduces fluctuations of the centroid 
parameter {zi}, it may be helpful in determining the matrix g of the AMD-QL 
model. Another possibility worth exploring is to incorporate the effects of two- 
particle correlations and their propagation 41 ) in a statistical manner. We hope that 
these challenges will help to deepen our microscopic understanding of fluctuation 
and dissipation phenomena in nuclear dynamics. 
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